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Abstract 

We establish a boundary Harnack principle for a large class of subordinate Brownian mo- 
tion, including mixtures of symmetric stable processes, in bounded K-fat open set (disconnected 
analogue of John domains). As an application of the boundary Harnack principle, we identify 
the Martin boundary and the minimal Martin boundary of bounded K-fat open sets with respect 
to these processes with their Euclidean boundary. 
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1 Introduction 



The boundary Harnack principle for nonnegative classical harmonic functions is a very deep result 
in potential theory and has very important applications in probability and potential theory. 

In [5] Bogdan showed that the boundary Harnack principle is valid in bounded Lipschitz domains 
for nonnegative harmonic functions of rotationally invariant stable processes and then in [28] Song 
and Wu extended the boundary Harnack principle for rotationally invariant stable processes to 
bounded K-fat open set. Subsequently Bogdan-Stos-Sztonyk [7j and Sztonyk [30] extended the 
boundary Harnack principle to symmetric (not necessarily rotationally invariant) stable processes. 
In a recent paper [6], Bogdan, Kulczycki and Kwasnicki proved a version of the boundary Harnack 
inequality for nonnegative harmonic functions of rotationally invariant stable processes in arbitrary 
open sets. 

By using some perturbation methods, the boundary Harnack principle has been generalized 
to some classes of rotationally invariant Levy processes including relativistic stable processes and 
truncated stable processes. These processes can be regarded as perturbations of rotationally in- 
variant stable processes and their Green functions on bounded smooth domains are comparable to 
their counterparts for rotationally invariant stable processes (see [9], [22], [12], [15], [16] and [T7]). 
This comparison of Green functions played a crucial role in the arguments of [12] , [16] and |17j . 

In this paper, we will show that, under minimal conditions, the boundary Harnack principle is 
valid for subordinate Brownian motions with characteristic exponents of the form $(<!;) = |£| a ^(|£| 2 ) 
for some a € (0,2) and some positive function t which is slowly varying at oo. Examples of this 
class of subordinate Brownian motions include, among others, relativistic stable processes and 
mixtures of rotationally invariant stable processes. The Green functions of subordinate Brownian 
motions considered here behave like up^wup^ near the origin. So these subordinate Brownian 
motions can not be regarded as perturbations of rotationally invariant stable processes in general 
and their Green functions in bounded smooth domains are not comparable to their counterparts 
for rotationally invariant stable processes. 

Our proof of the boundary Harnack principle will be similar to the arguments in [J] and [28] for 
rotationally invariant stable processes. One of the key ingredients is a sharp upper bound for the 
expected exit time from a ball which, in the case of stable processes, follows easily from the explicit 
formula for the Green function of a ball. However, in the present case, the desired upper bound 
is pretty difficult to establish. We rely on the fluctuation theory for real-valued Levy processes to 
accomplish this. 

The organization of this paper is as follows. In Section 2 we use the fluctuation theory for real- 
valued Levy processes to establish a nice upper bound on the expected exit time from an interval 
for a one-dimensional subordinate Brownian motion. In Section 3, we use the results of Section 2 
to establish the desired upper bound on the expected exit time from a ball for a multidimensional 
subordinate Brownian motion and an upper bound on the Poisson kernel of a ball. The proof of 
the boundary Harnack principle is given in Section 4 and in the last section we apply our boundary 
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Harnack principle to study the Martin boundary with respect to subordinate Brownian motions. 

In this paper we will use the following convention: the values of the constants r\^r%, . . . will 
remain the same throughout this paper, while the values of the constants c\,c%,... or C, C\, C*2, . . . 
might change from one appearance to another. The dependence of the constants on the dimension, 
the index a and the slowly varying function will not be mentioned explicitly, while the dependence 
of the constants on other quantities will be expressed using c(-) with the arguments representing 
the quantities the constant depends on. In this paper, we use ":=" to denote a definition, which 
is read as "is defined to be". f{t) ~ g(t), t — > {fit) ~ g{t), t — ► oo, respectively) means 
lim t _> f{t)/g{t) = 1 (lim^oo = 1, respectively). 

2 Some Results on One-dimensional Subordinate Brownian Mo- 
tion 

Suppose that W = (Wt ■ t > 0) is a one-dimensional Brownian motion with 



E 



e ii{w t -w ) 



and S = {St : t > 0) is a subordinator (a non-negative increasing Levy process) independent of W 
and with Laplace exponent (ft, that is 



E 



e -xs t 



e-**W, Vi,A>0. 



A C°° function g : (0, oo) — > [0, oo) is called a Bernstein function if {—l) n D n g < for every positive 
integer n. A Bernstein function g can be written in the following form 

/•oo 

g(\) = a + b\+ (1 - e' xt )n{dt) 
Jo 

where a, b > and fi is a measure on (0, oo) with f^°{lAt)fi{dt) < oo. fi is called the Levy measure 
of g. It is well known that a function g is the Laplace exponent of a (non-killed) subordinator 
if and only if g is a Bernstein function with lim^o = 0. A Bernstein function g is called a 
complete Bernstein function if its Levy measure fi has a completely monotone density with respect 
to the Lebesgue measure. For details on examples and properties of complete Bernstein functions, 
one can see [13], [23] or [27]. One important property we are going to use in this paper is that / is 
complete Bernstein is equivalent to that A// (A) is complete Bernstein. Throughout this paper we 
will assume that (j) is a complete Bernstein function such that 

0(A) = A Q / 2 ^(A) (2.1) 

for some a € (0, 2) and some positive function I which is slowly varying at oo. For concepts and 
results related to the slowly varying function, we refer our readers to [3]. 
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Using Corollary 2.3 of [26] or Theorem 2.3 of [21] we know that the potential measure of S has 
a decreasing density u. By using the Tauberian theorem (Theorem 1.7.1 in [3]) and the monotone 
density theorem (Theorem 1.7.2 in [3]), one can easily check that 

Let /i(t) be the density of the Levy measure of (p. It follows from Proposition 2.23 of [27] that 

a lit- 1 ) . . 

"«> ~ 2T(l-a/2)f ^ - ' °- ( ^ 3> 

The subordinate Brownian motion X = {Xt : t > 0) defined by = Ws t is a symmetric Levy 
process with the characteristic exponent 

$(0) = (j)(e % ) = \9\ a l(8 2 ), V0€M. 

Let x be the Laplace exponent of the ladder height process of X. (For the definition of the ladder 
height process and its basic properties, we refer our readers to Chapter 6 of pQ.) Then it follows 
from Corollary 9.7 of [TO] that 

x( A) = exp^-y o -^»de)=ex^-J o V-L-^ » dB j , VA > 0. (2.4) 

Under our assumptions, we have the following result. 

Proposition 2.1 The Laplace exponent \ of the ladder height process of X is a special Bernstein 
function, i.e., A/%(A) is also a Bernstein function. 

Proof. Define = X/4>(\). Let T be a subordinator independent of W and with Laplace 

exponent ip and let Y = (Yt : t > 0) be the subordinate Brownian motion defined by Yt = Wx t ■ Let 
^ be the characteristic exponent of Y. Then 

$(0)*(0) = cj){e 2 )^{e 2 ) = e 2 , V0 e m. 

Let p be the Laplace exponent of the ladder height process of Y . Then by (12.41) we have 

x(A W A) _ exp (- r wwyw . 



7T 



1 + 2 



= exp l-7o 1T02 d * 

/i p io g (0 2 A 2 ) \ 

Thus x is a special Bernstein function. □ 
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Proposition 2.2 // there are M > 1, 8 G (0, 1) and a nonnegative integrable function f on (0, 5) 
such that 

'i{\ 2 e 2 ) 



log 



then 



Proof. Using the identity 
we get easily from (I2.4|) that 



1 roo 

exp | - 

70 



</(#), V(fl,A)€(0,<5)x(M,oo), 

lim — = 1 

A^oc A«/ 2 (^(A 2 )) 1 /2 

i og ((W) 



(2.5) 
(2.6) 



1 + 



dd 



VA,/3 > 0, 



d6> 



£(A : 



2/)2^ 



1 



r/0 



By Potter's Theorem (Theorem 1.5.6 (1) in [3]), there exists Ao > 1 such that 
7(A 2 fl 2 ) 



log 



^(A 2 ) 



1 + 



1 log<9 
< 2- 



l + e 2 ' 



V(0,A) G I 1 , 00 ) x [Ao,oo). 



Thus by using the dominated convergence theorem in the first integral below, the uniform conver- 
gence theorem (Theorem 1.2.1 in [3j) in the second integral, and the assumption (|2.5p in the third 
integral, we have 



f°°, (i(\ 2 e 2 ) 

hm / log 
/o 



A— >oo 



^(a 2 ) ; i + 



d6 = lim 

A^oo 



oo pX 
+ / + 



log 



1 



£(A 2 ) J i + e 2 



de = o. 



□ 



In the case </>(A) = A a / 2 for some a £ (0, 2), the assumption of the proposition above is trivially 
satisfied. Now we give some other examples. 

Example 2.3 Suppose that a € (0,2) and define 

<£(A) = (A + l) a / 2 -l. 
Then </> is a complete Bernstein function which can be written as <j>(\) = X a ^ 2 £(X) with 

(A + If/ 2 - 1 



£(X) 



A°/ 2 



Using elementary analysis one can easily check that there is a nonnegative integrable function / on 
(0, 1) such that (|2.5jl is satisfied. 
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Example 2.4 Suppose < /3 < a < 2 and define 

0(A) = A a / 2 + A' 3 / 2 . 

Then <f) is a complete Bernstein function which can be written as 0(A) = A a/2 ^(A) with 

£( A) = l + A (/3-)/2. 

Using elementary analysis one can easily check that there is a nonnegative integrable function / on 
(0, 1) such that (|2.5p is satisfied. 

Example 2.5 Suppose that a £ (0, 2) and (5 € (0, 2 — a). Define 

0(A) = A a / 2 (log(l + A))' 3 / 2 . 

By using the facts that A and log(l + A) are complete Bernstein functions and properties of complete 
Bernstein functions (see |27J), one can easily check that <p is a complete Bernstein function. </> can 
be written as 0(A) = A a/2 ^(A) with 

^(A) = (log(l + A))^ 2 . 

To check that there is a nonnegative integrable function / on (0, 1) such that (|2.5p is satisfied, we 
only need to bound the function 

/ log(l + A 2 fl 2 y 
°H log(l + A2) 

for large A and small 9. We will consider two cases separately. Fix an M > 1 and a 9 < 1. 

(1) A > M, 9 < 1 and A > 1/9. In this case, by using the fact that for any a > the function 
x l— *■ jz^ is decreasing on (o, oo), we get that 

'log(l + A 2 6> 2 ) 



'° g V log(l + A' 2 ) ) ~ l0g 

< log 

< log 
= log 

(2) A > M, 6 < 1 and A < 1/9. In this case we have 

dog(l + A 2 # 2 ) 



Ml + A 2 ) 
log(l + \ 2 9 2 

log(l + A 2 ) 
log(# 2 ) + log(l + A 2 ) 

log(l + 9~ 2 ) 
\og{9 2 ) + log(l + 9- 2 
log(l + 9 2 ) - log(6> 2 ) 

log(l + 9 2 ) 



log 



log(l + A 2 ) 



log 



< 



< 



Ml + A 2 ) 
log(l + A 2 6< 2 

log(l + A 2 ) 
log(l + M 2 9 2 ) 

log(l + 9- 2 ) 
log(l + M 2 9 2 ) 
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Combining the results above one can easily check that there is a nonnegative integrable function / 
on (0, 1) such that (|2.5p is satisfied. 

Example 2.6 Suppose that a G (0,2) and (5 G (0, a). Define 

0(A) = A a / 2 (log(l + A))-^ 2 . 

By using the facts that A and log(l + A) are complete Bernstein functions and properties of complete 
Bernstein functions (see [27]), one can easily check that ^ is a complete Bernstein function. (f> can 
be written as 0(A) = A a/2 ^(A) with 

£(A) = (log(l + A))-^ 2 . 

Similarly to the example above, once can use elementary analysis to check that there is a nonnegative 
integrable function / on (0, 1) such that (|2.5p is satisfied. 

The method used to construct the complete Bernstein functions can be used to construct a whole 
class of complete Bernstein functions satisfying the assumptions of this paper. For instance, one 
can check that, for a G (0,2), (3 G (0,2 -a), functions like A a / 2 (log(l + log(l + \))f /2 } A a / 2 (log(l + 
log(l+log(l+A))))^ //2 , . . . are complete Bernstein functions satisfying the assumptions of this paper. 
Similarly, for any a G (0,2), (3 G (0,a), functions like A a / 2 (log(l + log(l + A)))~^/ 2 , A Q / 2 (log(l + 
log(l + log(l + A)))) -/3//2 , . . . are complete Bernstein functions satisfying the assumptions of this 
paper. 

In the remainder of this section we will always assume that the assumption of Proposition 12.21 is 
satisfied. It follows from Propositions 12. l l and 12.21 above and Corollary 2.3 of [26] that the potential 
measure V of the ladder height process of X has a decreasing density v. Since X is symmetric, we 
know that the potential measure V of the dual ladder height process is equal to V. 

In light of Proposition 12.21 one can easily apply the Tauberian theorem (Theorem 1.7.1 in [3]) 
and the monotone density theorem (Theorem 1.7.2 in [3]) to get the following result. 

Proposition 2.7 As x — > 0, we have 

x a/2 

~ r(l + a/2)(£(x- 2 )) 1 /2' 

x a/2-l 

V ^ ~ r(a/2)(^(x- 2 )) 1 /2- 
Proof. We omit the details. □ 

It follows from Proposition [2?2] above and Lemma 7.10 of [19] that the process X does not creep 
upwards. Since X is symmetric, we know that X also does not creep downwards. Thus if, for any 
a G R 3 we define 

T a = inf{t > : X t < a}, a a = inf{t > : X t < a}, 
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then we have 

P*(T a = O = 1, x > a. (2.7) 

Let G^ ' 00 ) {x, y) be the Green function of X^' 00 ^ , the process obtained by killing X upon exiting 
from (0, oo). Then we have the following result. 

Proposition 2.8 For any x, y > we have 

G^(x,y) = { + - < V, 

1 J x _yV{z)v{y + z - x)dz, x>y. 

Proof. By using (|2,7p above and Theorem 20 on page 176 of [I] we get that for any nonnegative 
function on / on (0, oo), 



E, 



(0,oo)\ 



oft 



JO 



v(z)f(x + z - y)v(y)dzdy . 



(2. 



where k is the constant depending on the normalization of the local time of the process X reflected 
at its supremum. We choose k = 1. Then 



E 3 



v ( z ) / v(y)f(x + y - z)dydz 



v (y) v(z)f(x + y - z)dzdy 



v(z) / v(w + z — x)f(w)dwdz 



px pa 

/ fW / 

JO Jx- 



f(w) / v(z)v(w + z — x)dzdw + I f(w) / v{z)v{w + z — x)dzdw . (2.9) 



On the other hand, 



E* 



(0,oo)x 



dt 



poo 
JO 



G^' Qo) (x,w)f(w)dw+ / G^{x,w)f{w)dw 



By comparing (|2.9p and (|2.10p we arrive at our desired conclusion. 



(2.10) 

□ 



For any r > 0, let G^ '^ be the Green function of X^°' r \ the process obtained by killing X upon 
exiting from (0, r). Then we have the following result. 

Proposition 2.9 For any R > 0, i/iere exists C = C{R) > swc/i that 

r a/2 £ a / 2 



G^ r \x,y)dy<C 



o 



(£( r -2))l/2 (£( x -2))l/2 



a € (0,r), r € (0, J2). 
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Proof. For any x € (0, r), we have 

G^ r \x,y)dy < f G^\x,y)dy 
n Jo 

x rx pr fx 

\ v(z)v(y + z — x)dzdy + / v(z)v(y + z — x)dzdy 
Jx-y Jx JO 



v{z) I v(y + z - x)dydz + / v{z) I v(y + z - x)dydz < 2V((0,r))V((0,x)). 

Jx-z JO Jx 

Now the desired conclusion follows easily from Proposition 12,71 and the continuity of V((Q,x)) and 

WVfi^- 2 )) 1 / 2 . 

□ 

As a consequence of the result above, we immediately get the following. 
Proposition 2.10 For any R > 0, there exists C = C(R) > such that 

r a/2 ( all ( _ T \a/2 \ 

3 Key estimates on Multi-dimensional Subordinate Brownian Mo- 
tions 

In the remainder of this paper we will always assume that d > 2 and that a G (0, 2). From now on 
we will assume that B = (Bt : t > 0) is a Brownian motion on M. d with 



E 



e i£-(B t -B ) 



e~ m \ V£eR d ,t>0. 



Suppose that S = (5* : t > 0) is a subordinator independent of B and that its Laplace exponent 
is a complete Bernstein function satisfying all the assumption of the previous section. More 
precisely we assume that there is a positive function I on (0, oo) which is slowly varying at oo such 
that 4>(X) = X a ^ 2 £(X) for all A > and that there is a nonnegative integrable function / on (0, 5) 
for some 5 > such that (|2.5j) holds. As in the previous section, we will use u(t) and /i(t) to denote 
the potential density and Levy density of S respectively. 

In the sequel, we will use X = (Xt : t > 0) to denote the subordinate Brownian motion defined 
by Xt = Bg t - Then it is easy to check that when d > 3 the process X is transient. In the case of 
d = 2, we will always assume the following: 

Al. The potential density u of S satisfies the following assumption: 

u(t)~cF- x t t^oo (3.1) 

for some constants c > and 7 < 1 . 
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Under this assumption the process X is also transient for d = 2. 

We will use G(x, y) = G(x — y) to denote the Green function of X. The Green function G of 
X is given by the following formula 

/>oo 

G{x) = / {A7rty d/2 e-W 2/( - u) u(t)dt, x € R d . 





Using this formula, we can easily see that G is radially decreasing and continuous in R d \ {0}. 

In order to get the asymptotic behavior of G near the origin, we need some additional assumption 
on the slowly varying function i. For any y,t,£ > 0, define 

K H {y,t) := \ TOP y < 
' CV ^ \ 0, y>\. 

We will always assume that 

A2. There is a £ > such that 

A u (y,t) < g(t), Vy,t>0, 

for some positive function g on (0, oo) with 

poo 

/ t (d-oi)/2-l e -t g ( t)dt < ^ 

Jo 

It is easy to check (see the proofs of Theorem 3.6 and Theorem 3.11 in [27]) that for the 
subordinators corresponding to Examples I2.3t42~6l Al and A2 are satisfied. 
Under these assumptions we have the following. 

Theorem 3.1 The Green function G of X satisfies the following 

aT({d-a)/2) 1 
[ ' ~ 2°+ 1 7r d / 2 r(l + a/2) \x\ d - a Z(\x\- 2 ) ' 11 ^ 

Proof. This follows easily from A1-A2, (12. 2j) above and Lemma 3.3 of [27]. We omit the details. 

□ 

Let J be the jumping function of X, then 

/•oo 

J(x)= / (4irt)- d / 2 e-^ 2/ ( 4t) n(t)dt, xGR d . 
Jo 

Thus J(x) = j(\x\) with 

roo 

j(r)= / (4irt)- d / 2 e~ r2 /^fi(t)dt, r > 0. 



o 
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It is easy to see that j is continuous in (0, oo). Since t \— > fj,(t) is decreasing, the function r \— > j(r) 
is decreasing on (0, oo). In order to get the asymptotic behavior of j near the origin, we need some 
additional assumption on the slowly varying function I. For any y, t, £ > 0, define 



TOT' y < f 



0, y > * 



We will always assume that 

A3. There is a £ > such that 

T u (y,t)<h(t), Vy,t>0 

for some positive function h on (0, oo) with 

t (d+a)/2-l e -t h ^ dt < ^ 



It is easy to check (see the proofs of Theorem 3.6 and Theorem 3.11 in [27]) that for the 
subordinators corresponding to Examples I2.3H2.61 A3 is satisfied. 

Theorem 3.2 The function j satisfies the following 



aT((d + a)/2) l(r~ 2 ) 
2i-«vr d / 2 r(l - a/2) ~ d 



J\ r ) ~ - i jt^ — ; rr — Ti — > r ^ 0. 



Proof. This follows easily from Al, A3, (|2.3|) above and Lemma 3.3 of [57]. We omit the details. 
□ 

For any open set D, we use td to denote the first exit time from D, i.e., td = inf{i > : 
X t i D}. Given an open set D C R d , we define XP(cu) = X t (u) if t < t d (uj) and X t D (uj) = d if 
t > td(w), where 9 is a cemetery state. We now recall the definition of harmonic functions with 
respect to X. 

Definition 3.3 Let D be an open subset ofW 1 . A function u defined on M. d is said to be 

(1) harmonic in D with respect to X if 

~E X [|u(X TS )|] < oo and u{x) = K x [u(X TB )] , x € B, 
for every open set B whose closure is a compact subset of D; 

(2) regular harmonic in D with respect to X if it is harmonic in D with respect to X and for each 
x £ D, 

u(x)=E x [u{X TD )}; 
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(3) harmonic for X if it is harmonic for X in D and vanishes outside D. 

In order for a scale invariant Harnack inequality to hold, we need to assume some additional 
conditions on the Levy density \x of S. We will always assume that 

A4. The Levy density fi of S satisfies the following conditions: there exists C\ > such that 

H(t) < Cm(t + 1), Vt>l. 

It follows from (12. 3p that for any M > there exists C2 > such that 

Hit) < C 2 //(2t), Vt€(0,M). 
Using A4, the display above and repeating the proof of Lemma 4.2 of [2T] we get that 

(1) For any M > 0, there exists C3 > such that 

j(r)<C 3 j(2r), Vre(0,M). (3.2) 

(2) There exists C4 > such that 

< C 4 j(r + 1), Vr>l. 

It is easy to check (see \27\) that for the subordinators corresponding to Examples I2.3H2.61 
A4 is satisfied. Therefore by Theorem 4.14 of [27] (see also [21]) we have the following Harnack 
inequality. 

Theorem 3.4 (Harnack inequality) There exist r\ £ (0, 1) and C > such that for every 
r € (0, ri), every xq € M. d , and every nonnegative function u on M. d which is harmonic in B(xo,r) 
with respect to X, we have 

sup u(y) < C inf u(y). 

y£B(x ,r/2) yeB(x ,r/2) 

For any bounded open set D in M d , we will use Gr>{x,y) to denote the Green function of 
X D . Using the continuity and the radial decreasing property of G, we can easily check that Gd is 
continuous in (D x D)\ {(x, x) : x € D}. 

Proposition 3.5 For any R > 0, there exists C = C(R) > such that for every open subset D 
with diam(D) < R, 

G D (x, y) < G(x, y) <C ^ V(x, y) G D x D. (3.3) 

t{\x — y\ z )\x — y\ a a 



12 



Proof. The results of this proposition are immediate consequences of Theorem 13.11 and the conti- 
nuity and positivity of l(r~ 2 )r d ~ a on (0, oo). □ 

The idea of the proof of the next lemma comes from [30J . 
Lemma 3.6 For any R > 0, there exists C = C(R) > such that for every r G (0, R) and xo G M d , 

r a/2 ( r _\ x _ £ |W 2 

Mtb M ] < C {l{r _ 2))l/2 (A(r _ |a _ a;o|) - 2))1/2 , * e B(x ,r). 
Proof. Without loss of generality, we may assume that x = 0. For x ^ 0, put Z t = Then 

\X\ 

Zt is a Levy process on M with 

E(e ieZt ) = E(e W W Xt ) = e ~*l e W 2 ), G R. 

Thus is the type of one-dimensional subordinate Brownian motion we studied in the previous 
section. It is easy to see that, if Xt G B(0,r), then \Zt\ < r, hence 

where f = inf{i > : \Zt\ > r}. Now the desired conclusion follows easily from Proposition 12.101 
□ 

Lemma 3.7 There exist r 2 G (0,ri] and C > sitc/i t/iai /or every positive r <r 2 and xq G M^, 



Proof. The conclusion of this Lemma follows easily from Theorem 13.21 above and Lemma 3.2 of 
[25]. □ 

Using the Levy system for X, we know that for every bounded open subset D and every / > 
and x € D, 

E x [f(X TD ); X TD _ + X TD ] = [_ f G D (x,z)J(z - y)dzf{y)dy. (3.4) 
For notational convenience, we define 

K D (x,y) := G D (x,z)J(z - y)dz, (x,y) G D x If. (3.5) 
Jd 

Thus (|3.4p can be simply written as 

E x [f(X TD ); X TD _ + X TD ] = [_ K D (x,y)f(y)dy. 

Jd c 

Using the continuities of Gd and J, one can easily check that Kr> is continuous on D x D°. 
As a consequence of Lemma 13 .6113. 71 and (|3,5p . we get the following proposition. 
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Proposition 3.8 There exist C\, C 2 > such that for every r G (0, r 2 ) and xq G R , 

r a/2 f r _ l x _ x lja/2 

K B{xo , r) (x,y) < C l3 (\y- X o\-r) ie{r _ 2))1/2 (£((r _ _ , q|) - 2))1/2 > (3-6) 



for all (x,y) G B(xQ,r) x B(x Q ,r) and 

r a c 

K B(x ,r){xo,y) > C 2 J(2(y - xq))-^-^, \fyeB(x ,r). (3.7) 



Proof. Without loss of generality, we assume xq = 0. For z G 5(0, r) and y G 5(0, r) 

|y| — r < |y| — |z| < |z — y\ < \z\ + \y\ < r + |y| < 2|y|. 
Thus by the monotonicity of J, 



J(2y) < J(z - y) < j(\y\ - r). (z, y) G 5(0, r) x B{0, r) . 
Applying the above inequality and Lemma l3.6H3.7l to (|3.5p . we have proved the proposition. □ 

Proposition 3.9 For every a G (0, 1) ; there exists C = C(a) > such that for every r G (0, r 2 ), 
xq G R d and x\,x 2 G B(xq, ar), 



K B(x ,r)(xi,y) < CK B ( Xo rf(x 2 ,y), y&B(x ,r) . 

Proof. This follows easily from the Harnack inequality (Theorem I3.4p and the continuity of 
Kb{x ,t)- For details, see the proof of Lemma 4.2 in [30]. □ 

As an immediate consequence of Theorem 13.21 we have 

Lemma 3.10 There exists r% G (0, r 2 ] such that for every y G M. d with \y\ < r%, 
aT((d + a)/2) £(\y\~ 2 ) < < 2 a aT((d + a)/2) £(\y\- 2 ) 



2 2 ~ a TT d / 2 T{\ - a/2) \y\ d+a ~ y , ~ ir d l 2 T(l - a/2) \y\ d+a ' 
The next inequalities will be used several times in the remainder of this paper. 

Lemma 3.11 There exist r 4 G (0, r%\ and C > such that 

g a/2 r a/2 

<C --, V0<s<r<4r 4 , (3. 



(£(s- 2 )) 1/2 : (e(r- 2 )) 1/r 

„l-a/2 r l-a/2 

TT^ < C m , V0 < s < r < 4r 4 , (3.9) 

(£(s- 2 )) 1/2 (£(r- 2 )) 1/2 

,1-a/a ^( S -2))V2 < Cr i-a/2 ^^-2^1/2^ V < s < r < 4r 4 , (3.10) 
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W01 ds , c Wgr. vo< r ^, ,,n) 

(£(s- 2 )) 1/2 U(r- 2 )) 1/2 , , 

{ S J ds ^ c [ r Ji > V0<r<4r 4 , (3.12) 

-j^ds <C^-^, V0<r<4r 4 , (3.13) 

£f s -2\ Mr-" 2 ) 

-±--±ds < C^—J-, V0<r<4r 4 (3.14) 



r 







and 



r fl a-l 



ds < C V0<r<4r 4 . (3.15) 



£(s- 2 ) " ^(r- 2 )' 



Proof. The first three inequalities follow easily from Theorem 1.5.3 of [3], while the last five from 
the 0- version of Theorem 1.5.11 of [3] ■ 1=1 

Proposition 3.12 For every a G (0, 1), there exists C = C{a) > such that for every r £ (0, r 4 ] 
and xq € 

a/2-d c^cc I _ I _ \-2\\l/2 

< C {e{r _ 2))l/2 ^| _ ^_ Vx e B(x 0l or), y 6 {r < |x -y| < 2r}. 



Proof. By Proposition 13.91 

K B(x (h r)(x,y) < ~2 / K B{x0tr) (w,y)dw 

' JB(xo,ar) 

for some constant ci = c\{a) > 0. Thus from Lemma 13.61 and (|3.6|) we have that 
^Bfso.O (^2/) - ~1 / / G B(x ,r)(w,z)J(z -y)dzdw 

r JB(x ,r) JB(x ,r) 



£2 



E *[TB(x ,r)K0 - 



c 3 r"/ 2 y (r-|z-x |) Q / 2 _ 

^ • ■ , -JB( m .rAl«r-\z-X \)->W AZ V) 



r d (£ (r -2 )) l/2 ( ^ ((r _ | Z _ X0 |)-2))l/2 ' 

for some constants C2 = 02(a) > and C3 = 03(0) > 0. Now applying Lemma f3. 101 we get 

c 4 W 2 - rf y (r-|z-xo|) a / 2 £(|z-y|- 2 ) 

&Bto,r)V,V) - {£{r -2 ) y/2j BMmr _ lz _ Xol) -2 )) l/2 \ z _ y]d+ a ** 

for some constant c 4 = c 4 (a) > 0. Since r — \z — xq\ < \y — z\ <3r< 3r 4 , from (|3.8|) we see that 

( r -\z-x \r/ 2 (\y-z\r/ 2 

< c 5 - 
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for some constant C5 > 0. Thus we have 



K B{xo>r) {x,y) < (£(r _ 2))1/2 



c 6 r a / 2 - d f (£(\z - y^ 2 )) 1 ' 2 



B(x ,r) \z-y\ d + a / 2 



-dz 



Cl r a / 2 - d r°° {e(s- 2 )) 1/2 

for some constants c§ = cq(o) > and cj = 07(a) > 0. Using (|3.1ip in the above equation, we 
conclude that 

K BM[ x,y) _ Wr _ 2))1/2 ( | y _ Xo |_ r)Q/2 
for some constant c$ = cs(a) > 0. □ 



4 Boundary Harnack Principle 

In this section, we give the proof of the boundary Harnack principle for X. 

Using an argument similar to the first part of the proof of Lemma 3.3 in [25] and using Lemma 
13.101 and (|3. 13j) - ()3. 14j) above we can easily get the following lemma. We skip the details. 

Lemma 4.1 There exists C > such that for any r £ (0, r^) and any open set D with D C B(0, r) 
we have 

P x (Xr D G-B(0,r) c ) < Cr- a i{r- 2 ) [ G D (x,y)dy, x E D n B(0, r/2). 

Jd 

Lemma 4.2 There exists C > such that for any open set D with B(A,nr) C D C B(0,r) for 
some r E (0, 7*4) and k £ (0, 1), we have that for every x E D\ B(A, nr), 

j D G D (x,y)dy < C^^- Q/2 ^^(l + ^|f^)^ {x TDmKr) eB(A,nr)). 

Proof. Fix a point x E D\ B(A,nr) and let B := B(A, ^p). Note that, by the harmonicity of 
G B (x, ■ ) in D \ {x} with respect to X, we have 

G D (x,A) > f K B (A,y)G D (x,y)dy > [ K B (A,y)G D (x,y)dy. 

JDnB c J Dr\B(A,^!f)c 

Since < \y — A\ < 2r for y E B(A, 2|l) c p jj an d j j s a decreasing function, it follows from 
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3.7p in Proposition 13.81 and Lemma 13.101 that 

/ Kr\a p 

G D (x,A) > Cl ^ / G D (x,y)J(2(y - A))dy 



' KT\—2\ 



^ C2K a r „ nKr ^-^ \ G D {x,y)dy, 

" ' 3kt\ 



for some positive constants c\ and C2. On the other hand, applying Theorem 13.41 we get 

/ G D (x,y)dy < c 3 / G D (x,A)dy < c 4 r d K d G D (x, A), 

Jb(a,^) Jb(a,^) 

for some positive constants C3 and C4. Combining these two estimates we get that 

J D G D (x,y)dy < c 5 ( rV + ^" Q ^Jpy) G D (x,A) (4.1) 

for some constant C5 > 0. 

Let O = D\B(A, Note that for any z G B(A, ^) and j/£fl, 2- 1 |y-z| < \y-A\ <2\y-z\. 
Thus we get from (f33j) that for 2 G 5(^4, 

c^K n {x,A) < K n {x,z) < c 6 K n {x,A) (4.2) 

for some cq > 1. Using the harmonicity of Gd(-,A) in D \ {A} with respect to X, we can split 
Gd(-,A) into two parts: 



G D (x,A) =E X [G D (X m ,A)} 

= E x [G D (X Tn ,A) : X Tn G B(A, ^)] + E, [g d pf Tfl , A) : X Tn G < \y - A\ < y} 
:=h + h. 



Using (|4.2p and (|3.3p . we have 

1 



/1 < c 6 K Q (x,A) [ G D {y,A)dy < c 7 K Q (x,A) [ 

JB(A,!f) Jb~ 



KA,f) JB{A^)\y-A\ d -i{\y-A\-^y 

for some constant c-j > 0. Since \y — A\ < 4r < 4r4, by (|3.8p . 

|y-^| Q/2 . (4r)"/ 2 

^ c 8 W9X (4.3) 



^(|y-^|~ 2 ) " °£((4r)" 2 ) 

for some constant cs > 0. Thus 

(4r) Q / 2 

I . A \rl — n/ r > , 



h < c 7 csKn(x,A) J ^^^^<^-^ )KaiXtA) 
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for some constant eg > 0. Now using (|4.2p again, we get 

1 



h < c 10 K a / 2 - d r a - d — / K n (x,z)dz, 



for some constant cio > 0. On the other hand, by (|3.3p . 

J 2 = / G D (y,A)P :E (X Tn Gdy) 

J { >f<\y- A \<mL } \y - A\ a a £(\y-A\ *) 
for some constant c\\ > 0. Using (|4.3|) . the above is less than or equal to 

for some constant c\2 > 0. Therefore 

G D (x,A) < c ^^ /2 ' dra ' % {{ l ) -2 ) ^ { X m e 
for some constant C13 > 0. Combining the above with ([4,ip . we get 



1 (, , <(( 



y)dy < c 14 r a k^I 2 -— — oT 1 + " 2 ' ' ) ^ ( Xr mB , , __. G ^ 



1? 



£((4r)- 2 ) V £((4r)- 2 ) / V t o\s ( a,^) v » 2 



for some constant C14 > 0. It follows immediately that 



D 



1 (, , mr 2 ) 

□ 



Combining Lemmas 14.1114.21 and using the translation invariant property, we have the following 

Lemma 4.3 There exists c\ > such that for any open set D with B(A,nr) C D C B(Q,r) for 
some r G (0, r 4 ) and k G (0, 1), we have that for every x G D f] B(Q, |), 

F x (X TD G B(Q, rf) < oy K - d -^ 2 j^-^ (l + |||py) P. (^ DXfl(A , Kr) € «r)) . 
Let A(x, a, b) := {y £R d : a < \y - x\ < b}. 



Lemma 4.4 Let D be an open set and < 2r < r 4 . For every Q G R and any positive function 

10_ u 
6 '' 6 



it vanishing on D c (~l -B(Q, V), i/iere is a o~ £ (^§r, -4-r) such that for any M G (1, 00) and x G 



DnB(Q,§r), 



u i X T D n B (Q,a)^ X r D nB(Q,a) &A(Q,a,M) < C — / J (y)u(y)dy 



/or some constant C = C{M) > independent of Q and u. 
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Proof. Without loss of generality, we may assume that Q = 0. Note that by (|3,12p 

my\-°r 2 ) 1/2 (\y\-vr a/ My)dyda 



x fr JA(0,a,2r) 



< ci 



<C 2 



A(0,-£r,2r) WO 



l((\y\-°r 2 ) 1/2 (\y\-vr a/2 d<Ju(y)dy 
£(8- 2 )V 2 8- a / 2 d8 \ u{y)dy 



\y\-fr 



lOr 



mv\-^r 2 ) l/2 (\v\~) l - a/2 »(v)dy 



X(0,lS!i,2r) 6 6 



for some positive constants ci and C2- Using (|3.1U|) . we get that there is a constant C3 > such 
that 



'((Ivl - ^r) 1/2 (M - ^ L ) 1 - tt/2 «(»)4/ < c 3 



A(0,if:,2r) 



i(\y\- 2 ) 1/2 \y\ 1 - a/2 u(y)dy, 



which is less than or equal to 



r l-a/2 



C 4 - 



£(\y\- 2 )u(y)dy 



^((2r)-2)V2y A(0ii o, )2r) ' 

for some constant C4 > by (|3.9p . Thus, by taking C5 > 6C2C4, we can conclude that there is a 
a G (^P-r, ^r) such that 

/ *((|yi-a)-y/^M-a)- a /yj,)d0 < C5 - r """ 

JA(0,a,2r) 

Let i£flfl 5(0, |r). Note that, since X satisfies the hypothesis H in [29], by Theorem 1 in 
1 29 1 we have 



£(( 2r )-2)i/2 7a ( o,^» 



£(|y|- 2 )«(y)dl/. (4.4) 



T Dnfl(0,(r) /> T DnB(0,o-) 



u(X 



T Dnfl(0,ir)/' T DnB(0,cr) 



X r D nB(0,cr) G ^(O) ";^)) r Z?nB(0,(T) = T B(0,a) 



E x 
= E X 
= E X 

<E X u{X );X €A(0,a,M) 

'A(0,<x,M) 

In the first equality above we have used the fact that u vanishes on D c nB(0, a). Since a < 2r < r±, 
from (|3.6p in Proposition 13.81 Proposition 13. 12l and Lemma 13.101 we have 



K B(o,a)(x,y)u(y)dy. 



E n 



< 



(■(> 



< X T D nB( ,.))i X T D nB (Q ^) £ A{0,a,M) 

0»/Z-d my \_ ^-2^1/2 



< 



f a a 

JA(Q ; <T,2r) ilia' 



a 



,a/2 



A(0,a,M) 

u(y)dy 



K B (o,a){x,y)u(y)dy 



+ C6 



A{0,a,2r) (l(a- 2 )) 1/2 (I! 

g/a ((7 -|x|r/ 2 
~ a i JTi =2va72 7777 1 i\-2\M/2 M ^ dy 

A(0,2r,A0 (A " )) 7 ~ M) )) 7 
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for some constant cq > 0. For y € A(0, 2r, M), ^\y\ < |y| — c and cr — \x\ < a < 2r. Thus by f)3.8[) 
and the monotonicity of j, we have 

n , \ o" a / 2 (<r- |x|) Q / 2 ,W, - a 

j(|y|- CT )77r^TTT7^7777 < 



(%- 2 ))V2 (^(( CT _ | x |)-2))i/2 - '^12 y £((2r)- 2 ) 
for some constant c-j > 0. Thus by (13.21) . we get 



. W 2 (<r- |dW 2 ,. r a 



(^(a- 2 )) 1 ^ (^(( CT _ | x |)-2))i/2 - °^ iy| V((2r) 

for some constant cs = c%{M) > 0. On the other hand, by (|3.8p and (|4.4|) . there exist positive 
constants Cg and Cio such that 

a «/2-d (£((| y |_ a)- 2 )) 1 ^ 



A(0,a,2r) (£(a- 2 )) 1/2 (\y\ ~ <?) a l 2 



-u(y)dy 



(o r \a/2 r —a/2 r 

< c 9 r- d ^ — — / £(\y\- 2 )u(y)dy 

W(2r)- 2 )) 1 / 2 (^((2r)- 2 )) 1 / 2 A(o,^» U ' ' ^ 



< c io777f^ / r K\y\- 2 )\yr d - a u(y)dy, 



which is less than or equal to 



£((2r) 2 ) 1 / J 4(o,3^ I ar) 



for some constants en > by Lemma 13.101 Hence 



u (^n 8M );^n»M G A(0,a,M) < cia p( Z. 2 s / in J{y)u{y)dy 

J H(2r) z j 7^(0 ior M) 



M(0,^,M) 

for some constant C12 = ci2(M) > 0. □ 



Lemma 4.5 Lei Z) be an open set. Assume that B(A,Kr) C Df]B(Q,r) for some < r < 2r^ and 
k £ (0, g-]. Suppose that u > is regular harmonic in D n -B(Q, 2r) loit/j respect to X and u = in 
(D c (~l -B(Q, 2r)) U B(Q, M) c . If w is a regular harmonic function with respect to X in D n -B(Q, r) 
such that 

j u(x), xGB(Q,f) c Up c nB(Q,r)), 
[ ] ~ I 0, xGA(Q,r,f), 

> w(A) > C n(x), xefln 5(Q, ^r) 

/or some constant C = C{M) > 0. 
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Proof. Without loss of generality, we may assume Q = and x G -D (1 -B(0, |r). The left hand side 
inequality in the conclusion of the lemma is obvious, so we only need to prove the right hand side 
inequality Since u is regular harmonic in D n -6(0, 2r) with respect to X and u = on B(0, M) c , 
we know from Lemma 14,41 that there exists a G (— , ^=£) such that 

u(x) =E X u(X T ); X T e A(0,a,M) < c x - ' / J{y)u(y)dy 

for some constant ci = ci(JVf) > 0. On the other hand, by (|3.7jl in Proposition 13,81 we have that 



u>(A) = / K DnB(0jr) (A,y)u(y)dy > / K B , AKr) {A,y)u(y)dy 



> c 2 / J(2(A-y)) ^ " 

7A(0,4f,Af) *((«>") ) 



'A(0,4p,M) 

for some constant c 2 > 0. Note that |y — A\ < 2|y| in j4(0, ^,M). Hence by the monotonicity of j 
and (I3T21 . 

«>(A) > c 2 Kr " / j(4|y|My)rfy > c 3 Kr " / J(y)u(y)dy 



for some constant C3 = C3(M) > 0. Therefore 

for some constant C4 = C4(M) > 0. □ 



We recall the definition of K-fat set from [28J. 

Definition 4.6 Let k G (0, 1/2]. VKe say i/tai an open se£ D in M. d is K-fat if there exists R > 
snc/i that for each Q G dD and r G (0, I? n B(Q,r) contains a ball B(A r (Q), kt). The pair 
(R, k) is called the characteristics of the K-fat open set D. 

Note that all Lipschitz domain and all non-tangentially accessible domain (see [14] for the 
definition) are K-fat. Moreover, every John domain is K-fat (see Lemma 6.3 in |20j). The boundary 
of a K-fat open set can be highly nonrectifiable and, in general, no regularity of its boundary can 
be inferred. Bounded K-fat open set may be disconnected. 

Since I is slowly varying at 00, we get the Carleson's estimate from Lemma 14.51 



Corollary 4.7 Suppose that D is a bounded K-fat open set with the characteristics (R,k). There 
exists a constant Ri such that if r < R\, Q G dD, u> is regular harmonic in D n B(Q, 2r) with 
respect to X and u = in (D c n B(Q, 2r)) U B(Q, M) c , then 



3 

u(A) > u(x), x G D n B(Q, -r) 



for some constant C = C{M) > 0. 
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The next theorem is a boundary Harnack principle for bounded K-fat open set and it is the main 
result of this section. Maybe a word of caution is in order here. The boundary Harnack principle 
here is a little different from the ones proved in [3] and [25] in the sense that in the boundary 
Harnack principle below we require our harmonic functions to vanish on the complement of some 
large open ball, which is a little weaker than assuming that it vanishes on the whole complement 
of the open set. However, this will not affect our application later since we are mainly interested 
in the case when the harmonic functions are given by the Green functions. 

Theorem 4.8 Suppose that D is a bounded K-fat open set with the characteristics (R,k). There 
exists a constant r§ := r^(D,a,l) < r\ A R such that if 2r < r§ and Q € dD, then for any 
nonnegative functions u, v in M. d which are regular harmonic in Dn B(Q, 2r) with respect to X and 
vanish in D c n B(Q, 2r) U B(Q, M) c , we have 

! u(A r (Q)) u{x) u(A r (Q)) Ta nnR(n r \ 
for some constant C = C(D, M) > 1. 

Proof. Since I is slowly varying at oo, there exists a constant r§ := r^(D,a,l) < r± A R such that 
for every 2r < r§, 

ma Jj^L ^ 2r r 2 ) ^(f)~ 2 ) ^r 2 ) \ < 2 , 45) 

V^CM -2 )' ^((4r)" 2 )' £((4r)- 2 )' £((2r)- 2 ) J ~ { ' 

Fix 2r < r§ throughout this proof. Without loss of generality we may assume that Q = and 
u(A r (0)) = v(A r (0)). For simplicity, we will write ^4 r (0) as A in the remainder of this proof. Define 
ui and U2 to be regular harmonic functions inDfl B(0, r) with respect to X such that 



u\(x) 



u(x), r < \x\ < Tp, 



0, x e £(0,f ) C U (D c n5(0,r)) 

and 



u(x), x e B(0, f ) c U {D c n B(0, r)), 



2 

0, r<\x\<%, 



u 2 (x) = 

and note that u = u\ + u 2 . If D n {r < \y\ < ^} is empty, then u\ = and the inequality (|4,9 



below holds trivially. So we assume D n {r < \y\ < 4^} is not empty. Then by Lemma S3 

u{y) - ciK ' a imSj u{A) ' yeDn m y } ' 
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for some constant c\ = c\{M) > 0. For i6Dn -6(0, £), we have 
u\{x) = E, 



3t* 

<^ ns(0>r) ) : X TDnBl0tr) £Dn{r<\y\<-} 



< ( sup <!/) Fxk nBM efln{r<|!,|<^} 

\DD{r<|»|<|: } / V 

< ( sup u(y))F x (x TDnB{0r) eB(0,ry 
\Dn{r<\ y \<f } J v 

- ci K ~ a jm^j u{A) Fx K-c^ G 5(05 r) 

Now using Lemma 14.31 and (j4.5j) we have that for a; € D n -6(0, |), 

ui(x) (4.6) 



< c 3 u(A) P, (Xr (Dnfl(0 , rmB(A , ¥ , E y )) (4.7) 
for some positive constants C2 = c^M) and C3 = cs(M, n). Since 2r < r±, Theorem 13.41 implies that 

u{y) > c 4 u(A), y G fl(i,|) 
for some constant C4 > 0. Therefore for iGDfl -6(0, |) 

u(s) = E, [u(X T(DnB{0rmB{A ^)] > c±u{A)W x (x T{DnB{0r ^ B{A ^ G B(A, ^)) . (4.8) 

Using (|4.7p . the analogue of (|4.8p for w and the assumption that u(A) = we get that for 

x eD(lB(0, §), 

ui(x) < c 3 P. (*r (DnB(0 , p)Afl(Ai¥) G S(A y )) < c 5 (4.9) 
for some constant C5 = cs(M, k) > 0. Since u = 0on -6(0, M) c , we have that for x G On -6(0, r), 



n 2 (x) = / K DnB{0r) (x,z)u(z)dz 



A(QM,M) JDnB(0,r) 



G D nB(o,r){x,y)J(y - z)dyu(z)dz. 



Let 



s{x) := / G DnB ^ r) (x,y)dy. 

JDnB(0,r) 



Note that for every y G .6(0, r) and zGB(0,| 



DnB(0,r) 

ir\c 
2 I ' 



— J jzr I < \z\ — r < |z| — |y| < \y — z\ < \y\ + \z\ < r + |z| < 2|z|. 
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So by the monotonicity of j, for every y G B(0,r) and z G -B(0, 

j(12|z|) < j(2kl) < J(y-z) < j(~\z\) < i(^N). 
Using (|3,2|) . we have that, for every y G -B(0, r) and 2 G A(0, M), 

^'(M) < -J(y-*) < c6i(N) 

for some constant C6 = C6(M) > 0. Thus we have 

, i^/iM < C7 , (4.10) 

for some constant C7 = ct(M) > 1. Applying (|4.10p to u and u and Lemma I3~5l to v and v%, we 
obtain for x G D n B(0, §), 

«(x) < or MA) ^ < c?^«(x) < ,«-JLJ_i-±J wW = ,«--JLJ_2 w(x)l 

(4.11) 

for some constant cs = c%(M) > 0. Combining (|4.9p and (|4.1ip and applying (|4.5p . we have 

u(sc) < cgv(x), ieflnB(o,p, 
for some constant eg = cg(M, k) > 0. □ 



5 Martin Boundary and Martin Representation 

In this section we will always assume that D is a bounded K-fat open set in W 1 with the character- 
istics (R, k). We are going to apply Theorem 14.81 to study the Martin boundary of D with respect 
to X. 

We recall from Definition ^. 6l that for each Q G dD and r G (0, R), A r (Q) is a point in Df]B(Q, r) 
satisfying B(A r (Q), nr) C DC\B(Q,r). From Theorem l4.81 we get the following boundary Harnack 
principle for the Green function of X which will play an important role in this section. Recall that 
r 5 < R is the constant defined in Theorem 14.81 



Theorem 5.1 There exists a constant c = c(D,a,£) > 1 such that for any Q G dD, r G (0, r§) 
and z,w G D \ B(Q, 2r), we have 

x G D (z,A r (Q)) G D (z,x) G D (z,MQ)) ^ D n B ( O — 
G D (w,A r (Q)) ~ G D {w,x) ~ G D (w,A r (Q)) ' V*' 2 

Since I is slowly varying at 00, there exists a positive constant r§ := Tq(k, I) < r$ such that for 
every 2r < re, 

2 - min ( -^r ! - max [~i(F*r' - ( } 
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Lemma 5.2 There exist positive constants c = c(D,a) and 7 = j(D, a) < a such that for any 
Q £ dD and r € (0,re), and nonnegative function u which is harmonic with respect to X in 
Dr\B(Q,r) we have 



u(A r (Q)) < c (-J U ^ » u(A {K/2)kr (Q)), k = 0,l,.... (5.2) 

Proof. Without loss of generality, we may assume Q = 0. Fix r < r& and let 

(k \ k 
-j r, A k := A Vh (0) and B k := B(A k ,r/ k+1 ), k = 0,l,.... 

Note that the B k s are disjoint. So by the harmonicity of u, we have 

k-l fc-i . 

u{A k ) > ^E Afc [u(Y TBk ) : Y TBk £B,j = K Bk (A k , z)u(z)dz. 

1=0 1=0 

Theorem 13.41 implies that 



/ K Bk (A k ,z)u(z)dz > cqu(Ai) / K Bk (A k ,z)dz 
Jb 1 Jb 1 

for some constant Co = co(d, a) > 0. Since dist(A k , Bi) < 2r]i, by (|3.7p in Proposition 13.81 and the 
monotonicity of j we have 

K Bk (A k ,z) > d J(2(A fc - z)) > Cl J(4tb) (,7fc+l) ° z € B,. 

mVk+i) 2 ) ) 

Applying Lemma 13.101 and (|5.ip , for every z in Z?/ we get 

K Bk (A k ,z) > c 2 7— . .., ... ^2T777Z T^oT ^ 2c 2 o 



£((^ +l) -2) i(( %+1 )- 2 ) - v 8 ; ^((% + i)- 2 ) 

for some constant C2 = C2(<i, a, £) > 0. Thus we have 

/ & B {A k , z)dz > c 3 777777 r~9T' ze^l 

7b ; ^((%+i)" 2 ) 

for some constant C3 = C3(d, a,£) > 0. Therefore, 

fc-i 

( Vk y a u (A k )£({ Vk+1 r 2 ) > c^imY^uiA^i^y 2 ) 



1=0 

■vfc-1 



for some constant C4 = C4(d, a, k,£) > 0. Let := (rj k ) a u(A k )£( (^^2 ) so that > C4 Xw=o 
By induction, one can easily check that a k > 05(1 + Ci/2) k ao for some constant C5 = 05(^,0;) > 0. 
Thus, with 7 = a - ln(l + ^){hx{2/ n))' 1 , we get 

{2\ lk I ((K/2)~ 2 ( fc+1 V~ 2 )) 
U(A(Q)) < C (^-J l(( K / 2 )-2 r ~ 2 )) U ( A (K/V k r(Q))- 

Applying (|5.ip . we conclude that (|5.2p is true. □ 
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Lemma 5.3 Suppose Q E dD and r G (0, r^). If w E D\B(Q,r), then 

G D (A r (Q),w) > c ^ [ J{\{z-Q))G D (z,w)dz 

for some constant c = c(D, a, t) > 0. 

Proof. Without loss of generality, we may assume Q = 0. Fix w £ D \ B(0, r) and let A := ^4 r (0) 
and u(-) := Gd(-,w). Since u is regular harmonic in D n -6(0, (1 — n/2)r) with respect to X, we 
have 



> E A [u [X TDnBmi _ K/2)r) j ;X TDnB{0 (1 _ K/2)r) G B(0, r) c 
K D nB(o,(i- K /2)r)iA z)u(z)dz 

B(0,r) c 

/ G DnB (p,(i- K /2)r) (-4, y)J(y~ z)dyu(z)dz. 

>B(0,r) c JDnB(p,(l-K,/2)r) 

Since B(A, nr/2) C D n -6(0, (1 — k/2)t), by the monotonicity of the Green functions, 

GDnB(0,(l- K /2)r)(Ay) > G f B(A,/cr/2)(Ay)) y G B(A, ftr/2). 

Thus 



> / / G B{A)Kr / 2 ){A,y)J(y - z)dyu{z)dz 

JB(0,r) c J B(A,nr/2) 

K B (A, K r/2)(A Z)u(z)dz, 

B(0,r) c 

which is greater than or equal to 

{Kr/2) c 

Ib(o,t)°" ""l{{nr/2) 

for some positive constant c% = c\(d,a,l) by (|3.7[) in Proposition 13.81 Note that \z — A\ < 2\z\ for 
z G B(0,r) c . Let M :=diam(D). Hence 

u(A)>c 2 -- — 77^-27/ u(z)J(4:z)dz > c 3 — — -—^ / u(z)J(-z)dz (5.3) 

£({Kr/2) 2 ) J A (o,r,M) £{{Kr/2) 2 ) Ja^m) 2 

for some constant C3 = C3(d, a, M) > 0. We have used (|3.2p in the last inequality above. □ 



ci j J(2(z - A^ ^ JJ u{z)dz 



Lemma 5.4 There exist positive constants c\ = ci(D,a,l) and C2 = C2(D,a,l) < 1 such that for 
any Q G dD, r G (0, r§) and w G D\ B(Q, 2r j k), we have 



E.r 



G D (X T , w) : X T G B(Q, r) c < ci Cg Gd(x, w), x (£ D n B k 



where B k := B(Q, (K/2) k r), k = 0, 1, . . . . 
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Proof. Without loss of generality, we may assume Q = 0. Fix r < r§ and w G D \ B(0,4r). Let 
rj k := (K/2) k r, B k := B(0,r] k ) and 



u k {x) := E x G D (X TDnB w);X TDnB £ B(0,r 



, xeDnB k . 



Note that for x £ D D B k+1 



u k+1 (x) 



= IE. 

= E x 

= E a . 

< E T 



G D (X Tj3nBfe+i ,u;); TDnB fc+1 = T Dri B k , X TDnBk+i G 5(0^)' 

Gfl^nB^w); T DnB k+1 = T D nB k , X r DnBk G 5(0^)' 

Go^n,,™); X TDnBh G -6(0, r) 



Thus 



Let A k := ^4^(0) and M :=diam(D). Since Gd( • ,w) is zero on D c , we have 
u k (A k ) = E Ak \G D (X TDnB w); X TDnBh G A(0, r, M) 



(5.4) 



< E 



G_d (X Tg , w) ; X T „ 6 4(0, r,M) 



< 



'A(0,r,M) 

Since r < r4, by (|3.6p in Proposition 13.81 we get that for z G A(0, r, M), 

fa* - |>U r 12 



K Bk (A k ,z)G D (z,w)dz. 



K Bk (A k ,z) < ci J(\z\ - rj k ) 



(^- 2 ))V2 (£(( m - |A fc |)-2))l/2 



for some constant ci = ci(D, a) > and fc = 1, 2, . . . . Since r/^ — \A k \ < rj k < r§, from ()3.8p we see 
that 

(m-\A k \) a / 2 



< c- 



a/2 



Thus 



K Bk (A k ,z) < c 2 J(\z\ - r] k ) 



V k 



?(v k 2 ) 



for some constant c% = C2(D, a,£) > and k = 1, 2, ... . Therefore by the monotonicity of j 



u k (A k ) < c 2 



V k 



£(v k 2 ) Ja(o,t,m) 



From Lemma 15.31 we have 

G D (A ,w) > c 3 



J(-z)G]j{z,w)dz, /c = 1,2, 



J(-z)Gd(z, w)dz 



£((W2)- 2 ) J A (0,r,M) 2 



(5.5) 



(5.6) 
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for some constant C3 = cs(D,a,£) > 0. Therefore (|5.5p and (|5,6p imply that 

nto I ((n/2)~ 2 r- 2 ) 
l) i{{K/2)' 2k r- 2 ) 



for some constant C4 = c±{D,a,t) > 0. On the other hand, using Lemma 15.21 we get 
G D (A ,w) < c 5 (1Y H e -^-^pl GD (A k ,w) 



for some constant C5 = c^{D,a) > 0. Thus by (|5.1 



• 2\ -fc(«-7) 



for some constant C6 = cq(D, a) > and fe = 1, 2, . . . . By Theorem 15 .1\ we have 

< TT-, \ < c 6 „ 7 a \ < C 4 C 5 C 6 



G D (x,w) G D (x,w) G D (A k _i,w) \k / 

for k = 1,2,... . □ 

Let € -D be fixed and set 

M D (x,y) := ^ D ^ X,V \ , x,yeD,y^x . 

GD(xo,y) 

Md is called the Martin kernel of D with respect to X. 

Now the next theorem follows from Theorem 15.11 and Lemma 15.41 (instead of Lemma 13 and 
Lemma 14 in [1] respectively) in very much the same way as in the case of symmetric stable 
processes in Lemma 16 of [I] (with Green functions instead of harmonic functions). We omit the 
details. 

Theorem 5.5 There exist positive constants R±, M\, c and (3 depending on D, a and I such that 
for any Q € dD, r < R\ and z E D \ B(Q, Mir), we have 

\M D {z,x)-M D (z,y)\ < c ( \^lA \ x ,yeDnB(Q,r). 



In particular, the limit limD^ y ^ w Mq{x, y) exists for every w € dD. 

There is a compactification D M of D, unique up to a homeomorphism, such that Mz)(x,y) has 
a continuous extension toDx (D M \ {xq}) and Md(-, z\) = Md(-, Z2) if and only if Z\ = z%. (See, 
for instance, [15].) The set d M D = D M \D is called the Martin boundary of D. For z G d M D, set 
Md(-, z) to be zero in D c . 
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A positive harmonic function u for X D is minimal if, whenever v is a positive harmonic function 
for X D with v < u on D, one must have u = cv for some constant c. The set of points z G d M D 
such that Md(-, z) is minimal harmonic for X B is called the minimal Martin boundary of D. 

For each fixed z G 3D and x G D, let 

M D (x,z) := lim M D (x,y), 

which exists by Theorem 15,51 For each 2 G <9F>, set Mb(i, z) to be zero for x G -D c . 

Lemma 5.6 For every z G <9F> and B C B C D, Md(X tb , z) is ¥ x -integrable. 

Proof. Take a sequence {z m } m >i C D\B converging to 0. Since -/Vf£)(',z m ) is regular harmonic 
for X in by Fatou's lemma and Theorem 15.5^ 



E^ [M D (X rB ,z 



E, 



lim Mo (X Tg , z m ) < hm inf Md(x, z m ) = Mu(x,z) < 00. 



□ 



Lemma 5.7 For ewery z € and x £ D, 

M D (x, z) = E x M D (x? B(x ry z^ , for every < r < r 6 A -p D (x). 
Proof. Fix z G dD, x G D and r < r 6 A ^pd(^) < R- let 

Note that 



(5.7) 



-I r and z m := A„ ro (0), m = 0, 1,.... 



B(«m, Wi) C 5(2, -r/ m ) n D C B(z, r/ m ) n F> C B(z, r)flD C D\ B{x, r) 



for all m > 0. Thus by the harmonicity of Md(-, z m ), we have 

M D (x,z m ) = E,, [m d (x t 



TB(a:,r) ' ' 



On the other hand, by Theorem 15. 1\ there exist constants mo > and c\ > such that for 
every w € D\ B{z,rf m ) and y £ D (~) B(z,7] m+1 ), 



M D (w,z„ 



Letting y — » z G <9.D we get 



GD(w,z m ) Gr>(w,y) 

< ci— = ciM D [w,y), m>m . 



G D (x ,z, 



Gd(xo,v) 



M D (w,z m ) < ciM D (w,z), m>m , 



(5i 
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for every w £ D\ B(z, r) m ). 

To prove (|5.7p . it suffices to show that {Md(X Tb , r) , z m ) : m > mo} is P x -uniformly integrable. 
Since Mz>(X Tg , , , 2) is P x -integrable by Lemma ESI for any e > 0, there is an N > 1 such that 

Md (X T ,z);M D (x TB ^ r) , z) > N /a 



E* 

Note that by ([53]) and §x9§ 

E x [M D (x TB(x rV z m 

< a e, [m d (x TB(xr) , z) ■ Cl M D (x TB(xr) , z) > 

By (|3.6p in Proposition 13.81 we have for m > mo, 



< 



4ci 



(5.9) 



; M D \X TB{xr) ,z m ) > N and X TB(xr) £ D\B(z,m, 

No 



e e 
<C1 4^ = 4- 



E, 



< c 2 



„a/2 



DnB(z, Vm ) 

(r - H) a / 2 



M D (w,z m )j(|w -x| -r) iw 2 m/2 
DnB{ z , Vm ) Mr 2 )) l l 2 {£{{r - \w\) 2 ))V2 



M D (w, z m )K B ( x ^(x,w)dw 
dw 



for some c 2 = C2(d, a, /) > 0. Since \w — x\ > \x — z\ — \z — w\ > Pd( x ) — "Hm > 2r — r = r, usin^ 
the monotonicity of J and (|3,8p to the above equation, we see that 



E, 



< c 3 j(r) 

< C4 



M D (w,z m )dw 



£( r 2 ) JDnB(z, Vm ) 

M D (w,z m )dw = c i G D (xo,z rn y 1 



B(z,r) m ) 



(5.10) 



G D (w,z m )dw (5.11) 



'B(z,r)m) 

for some C3 = c^{D,a,£) > and C4 = c±(D, a, £,r) > 0. Note that, by Lemma 15.21 there exist 
C5 = c^(D, a, t, uiq) > 0, cq = cq(D, a, £, mo, r) > and 7 < a such that 



G D (x ,z m ) 1 < 



C5 I ^ 



< c 6 

On the other hand, by (|3.3p 



KX-7m £ (( K /2)- 2 ( m+1 )( K /2)- 2m °r- 2 )) 
^((K/2)- 2 (K/2)~ 2m o r -2) 

£ f(/€/2)- 2m (K/2)~ 2 ( m » +1 V- 2N ) . 



G D (x ,z n 



K\ 



(5.12) 



/ G D (w,z m )dw < c 7 

J B(z,r] m ) JE 



dw 



< c 8 



Ui: lll :* hn )t(\l»-Zm\ 2 )\w-Z m \ d a 
ds < Cg 



2rj m gCe—1 



(5.13) 



o £(s- 2 ) ~ s £((2r ?m )" 2 ) 

In the last inequality above, we have used (|3.15p . It follows from (|5.11j) - (|5.13j) that there exists 
cio = ao(D, a, £, mo, r) > such that 

Kx(a- 7 )m £ ((K/2)- 2m (K/2)- 2 ( mo+1 )r- 2 ) 



E, 



M D (X? B{xr) ,z m );X TB(xr) G Dfl B(z, 2r/m) 



< ao 



£{( K /2)- 2m (2r)- 2 ) 
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M d (x TB(xrV z m ) ; M D (x TB{xr) ,z m \ > V 



Since I is slowly varying at oo, we can take N = N(e, D, mo, r) large enough so that for m > N, 

IE. 

< E x \M D [X TB(x r) , z m J ; X TB(x r) £ D P\ B(z, 2r/m) 

+E a . [Md (X rs(:Cir) ,z m ) ; Mc (^,,4,) > N and X Tfl(a>ir) G -D \ B(z, 2r/m) 

/«x(a- 7 )m £((K/2)- 2m ( K /2)- 2 ( m o +1 )r- 2 ) e 

< 010 UJ £(( K /2)" 2m (2r)- 2 ) + 4 < £ ' 

As each Mr>(X Tg , . , z m ) is P x -integrable, we conclude that {Md(I Tb(j , , 2 m ) : m > mo} is uni- 
formly integrable under ¥ x . □ 



Using the fact that P x (X Tu G dU) = for every smooth open set U (Theorem 1 in [29]), one 
can follow the proof of Theorem 2.2 of [8] or the proof of Theorem 4.8 of [T7] and show that The 
two lemmas above imply that Md(-, z) is harmonic for X. We skip the details. 

Theorem 5.8 For every z G dD, the function x 1— > Md(-, z) is harmonic in D with respect to X . 

Recall that a point z G dD is said to be a regular boundary point for X if IP 2 (t£> = 0) = 1 and 
an irregular boundary point if P 2 (t£> = 0) = 0. It is well known that if z G <9D is regular for X, 
then for any x £ D, Gd(x, y) — ► as y — > z. 

Lemma 5.9 fi) If z,w G z ^ w and w is a regular boundary point for Y , then Mr>{x, z) — > 
OS x -> w. 

flJj T/ie mapping (x,z) 1— > Md(x,z) is continuous on D x <9D. 

Proof. Both of the assertions can be proved easily using our Theorems 15.11 and 15.51 We skip the 
proof since the argument is almost identical to the one on page 235 of [5] . □ 

Lemma 5.10 Suppose that h is a bounded singular a-harmonic function in a bounded open set D. 
If there is a set N of zero capacity such that for any z G dD \ N , 



lim h(x) = 0, 



then h is identically zero. 



Proof. Take an increasing sequence of open sets {D m } m >i satisfying D m C D m+ i and Um=l D m = 
D. Set r m = TD m . Then r m j tjj and lirrim^oo X Tm = X TD by the quasi-left continuity of X. Since 
N has zero capacity, we have 

P X (X TD G N) = 0, x G £>. 
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Therefore by the bounded convergence theorem we have for any x £ D, 



h{x) 



lim E x (h(X Tm ),T m < t d ) 

<i— >oo 

lim E x (h(X Tm )l dD \ N (X TD );T m < t d ) = 0. 



rn- 



rn- 



□ 



So far we have shown that the Martin boundary of D can be identified with a subset of the 
Euclidean boundary dD. 

If / is the set of irregular boundary points of D for X, then / is semi-polar by Proposition II. 3. 3 
in [2] j which is polar in our case (Theorem 4.1.2 in [H]). Thus Cap(I) = 0. Using this observation 
and the above lemma, now we can follow the proof of Theorem 4.1 in [2S| and show the following 
theorem, which is the main result of this section. 

Theorem 5.11 The Martin boundary and the minimal Martin boundary of D with respect to X 
can be identified with the Euclidean boundary of D. 

As a consequence of Theorem 15.111 we conclude that for every nonnegative harmonic function 
h for X D , there exists a unique finite measure \i on dD such that 



J 3D 

jj, is called the Martin measure of h. 
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